Application of the penalty coupling method for the analysis of blood vessels

Due to the significant health and economic impact of blood vessel diseases on modern society, its analysis is becoming of increasing importance for the medical sciences. The complexity of the vascular system, its dynamics and material characteristics all make it an ideal candidate for analysis through fluid structure interaction (FSI) simulations. FSI is a relatively new approach in numerical analysis and enables the multi-physical analysis of problems, yielding a higher accuracy of results than could be possible when using a single physics code to analyse the same category of problems. This paper introduces the concepts behind the Arbitrary Lagrangian Eulerian (ALE) formulation using the penalty coupling method. It moves on to present a validation case and compares it to available simulation results from the literature using a different FSI method. Results were found to correspond well to the comparison case as well as basic theory

Contrast-enhanced micro-CT imaging in murine carotid arteries : a new protocol for computing wall shear stress

Background: Wall shear stress (WSS) is involved in the pathophysiology of atherosclerosis. The correlation between WSS and atherosclerosis can be investigated over time using a WSS-manipulated atherosclerotic mouse model. To determine WSS in vivo, detailed 3D geometry of the vessel network is required. However, a protocol to reconstruct 3D murine vasculature using this animal model is lacking. In this project, we evaluated the adequacy of eXIA 160, a small animal contrast agent, for assessing murine vascular network on micro-CT. Also, a protocol was established for vessel geometry segmentation and WSS analysis. Methods: A tapering cast was placed around the right common carotid artery (RCCA) of ApoE(-/-) mice (n = 8). Contrast-enhanced micro-CT was performed using eXIA 160. An innovative local threshold-based segmentation procedure was implemented to reconstruct 3D geometry of the RCCA. The reconstructed RCCA was compared to the vessel geometry using a global threshold-based segmentation method. Computational fluid dynamics was applied to compute the velocity field and WSS distribution along the RCCA. Results: eXIA 160-enhanced micro-CT allowed clear visualization and assessment of the RCCA in all eight animals. No adverse biological effects were observed from the use of eXIA 160. Segmentation using local threshold values generated more accurate RCCA geometry than the global threshold-based approach. Mouse-specific velocity data and the RCCA geometry generated 3D WSS maps with high resolution, enabling quantitative analysis of WSS. In all animals, we observed low WSS upstream of the cast. Downstream of the cast, asymmetric WSS patterns were revealed with variation in size and location between animals. Conclusions: eXIA 160 provided good contrast to reconstruct 3D vessel geometry and determine WSS patterns in the RCCA of the atherosclerotic mouse model. We established a novel local threshold-based segmentation protocol for RCCA reconstruction and WSS computation. The observed differences between animals indicate the necessity to use mouse-specific data for WSS analysis. For our future work, our protocol makes it possible to study in vivo WSS longitudinally over a growing plaque

The definition of low wall shear stress and its effect on plaque progression estimation in human coronary arteries

Wall shear stress (WSS), the frictional force of the blood on the vessel wall, plays a crucial role in atherosclerotic plaque development. Low WSS has been associated with plaque growth, however previous research used different approaches to define low WSS to investigate its effect on plaque progression. In this study, we used four methodologies to allocate low, mid and high WSS in one dataset of human coronary arteries and investigated the predictive power of low WSS for plaque progression. Coronary reconstructions were based on multimodality imaging, using intravascular ultrasound and CT-imaging. Vessel-specific flow was measured using Doppler wire and computational fluid dynamics was performed to calculate WSS. The absolute WSS range varied greatly between the coronary arteries. On the population level, the established pattern of most plaque progression at low WSS was apparent in all methodologies defining the WSS categories. However, for the individual patient, when using measured flow to determine WSS, the absolute WSS values range so widely, that the use of absolute thresholds to determine low WSS was not appropriate to identify regions at high risk for plaque progression

Lipid-rich Plaques Detected by Near-infrared Spectroscopy Are More Frequently Exposed to High Shear Stress

High wall shear stress (WSS) and near-infrared spectroscopy (NIRS) detected lipid-rich plaque (LRP) are both known to be associated with plaque destabilization and future adverse cardiovascular events. However, knowledge of spatial co-localization of LRP and high WSS is lacking. This study investigated the co-localization of LRP based on NIRS and high WSS. Fifty-three patients presenting acute coronary syndrome underwent NIRS-intravascular-ultrasound (NIRS-IVUS) imaging of a non-culprit coronary artery. WSS was obtained using WSS profiling in 3D-reconstructions of the coronary arteries based on fusion of IVUS-segmented lumen and CT-derived 3D-centerline. Thirty-eight vessels were available for final analysis and divided into 0.5 mm/45° sectors. LRP sectors, as identified by NIRS, were more often colocalized with high WSS than sectors without LRP. Moreover, there was a dose-dependent relationship between lipid content and high WSS exposure. This study is a first step in understanding the evolution of LRPs to vulnerable plaques. [Figure not available: see fulltext.

Protein synthesis rates of muscle, tendon, ligament, cartilage, and bone tissue in vivo in humans

Skeletal muscle plasticity is reflected by a dynamic balance between protein synthesis and breakdown, with basal muscle tissue protein synthesis rates ranging between 0.02 and 0.09%/h. Though it is evident that other musculoskeletal tissues should also express some level of plasticity, data on protein synthesis rates of most of these tissues in vivo in humans is limited. Six otherwise healthy patients (62±3 y), scheduled to undergo unilateral total knee arthroplasty, were subjected to primed continuous intravenous infusions with L-[ring-13C6]-Phenylalanine throughout the surgical procedure. Tissue samples obtained during surgery included muscle, tendon, cruciate ligaments, cartilage, bone, menisci, fat, and synovium. Tissue-specific fractional protein synthesis rates (%/h) were assessed by measuring the incorporation of L-[ring-13C6]-Phenylalanine in tissue protein and were compared with muscle tissue protein synthesis rates using a paired t test. Tendon, bone, cartilage, Hoffa’s fat pad, anterior and posterior cruciate ligament, and menisci tissue protein synthesis rates averaged 0.06±0.01, 0.03±0.01, 0.04±0.01, 0.11±0.03, 0.07±0.02, 0.04±0.01, and 0.04±0.01%/h, respectively, and did not significantly differ from skeletal muscle protein synthesis rates (0.04±0.01%/h; P>0.05). Synovium derived protein (0.13±0.03%/h) and intercondylar notch bone tissue protein synthesis rates (0.03±0.01%/h) were respectively higher and lower compared to skeletal muscle protein synthesis rates (P<0.05 and P<0.01, respectively). Basal protein synthesis rates in various musculoskeletal tissues are within the same range of skeletal muscle protein synthesis rates, with fractional muscle, tendon, bone, cartilage, ligament, menisci, fat, and synovium protein synthesis rates ranging between 0.02 and 0.13% per hour in vivo in humans

Coronary fractional flow reserve measurements of a stenosed side branch: A computational study investigating the influence of the bifurcation angle

Background: Coronary hemodynamics and physiology specific for bifurcation lesions was not well understood. To investigate the influence of the bifurcation angle on the intracoronary hemodynamics of side branch (SB) lesions computational fluid dynamics simulations were performed. Methods: A parametric model representing a left anterior descending-first diagonal coronary bifurcation lesion was created according to the literature. Diameters obeyed fractal branching laws. Proximal and distal main branch (DMB) stenoses were both set at 60%. We varied the distal bifurcation angles (40°, 55°, and 70°), the flow splits to the DMB and SB (55%:45%, 65%:35%, and 75%:25%), and the SB stenoses (40, 60, and 80%), resulting in 27 simulations. Fractional flow reserve, defined as the ratio between the mean distal stenosis and mean aortic pressure during maximal hyperemia, was calculated for the DMB and SB (FFRSB) for all simulations. Results: The largest differences in FFRSB comparing the largest and smallest bifurcation angles were 0.02 (in cases with 40% SB stenosis, irrespective of the assumed flow split) and 0.05 (in cases with 60% SB stenosis, flow split 55%:45%). When the SB stenosis was 80%, the difference in FFRSB between the largest and smallest bifurcation angle was 0.33 (flow split 55%:45%). By describing the PSB-QSB relationship using a quadratic curve for cases with 80% SB stenosis, we found that the curve was steeper (i.e. higher flow resistance) when bifurcation angle increases (P=0.451*Q+0.010*Q 2 and P=0.687*Q+0.017*Q 2 for 40° and 70° bifurcation angle, respectively). Our analyses revealed complex hemodynamics in all cases with evident counter-rotating helical flow structures. Larger bifurcation angles resulted in more pronounced helical flow structures (i.e. higher helicity intensity), when 60 or 80% SB stenoses were present. A good correlation (R2=0.80) between the SB pressure drop and helicity intensity was also found. Conclusions: Our analyses showed that, in bifurcation lesions with 60% MB stenosis and 80% SB stenosis, SB pressure drop is higher for larger bifurcation angles suggesting higher flow resistance (i.e. curves describing the PSB-QSB relationship being steeper). When the SB stenosis is mild (40%) or moderate (60%), SB resistance is minimally influenced by the bifurcation angle, with differences not being clinically meaningful. Our findings also highlighted the complex interplay between anatomy, pressure drops, and blood flow helicity in bifurcations

2 days versus 5 days of postoperative antibiotics for complex appendicitis:a pragmatic, open-label, multicentre, non-inferiority randomised trial

Background: The appropriate duration of postoperative antibiotics for complex appendicitis is unclear. The increasing global threat of antimicrobial resistance warrants restrictive antibiotic use, which could also reduce side-effects, length of hospital stay, and costs. Methods: In this pragmatic, open-label, non-inferiority trial in 15 hospitals in the Netherlands, patients with complex appendicitis (aged ≥8 years) were randomly assigned (1:1) to receive 2 days or 5 days of intravenous antibiotics after appendicectomy. Randomisation was stratified by centre, and treating physicians and patients were not masked to treatment allocation. The primary endpoint was a composite endpoint of infectious complications and mortality within 90 days. The main outcome was the absolute risk difference (95% CI) in the primary endpoint, adjusted for age and severity of appendicitis, with a non-inferiority margin of 7·5%. Outcome assessment was based on electronic patient records and a telephone consultation 90 days after appendicectomy. Efficacy was analysed in the intention-to-treat and per-protocol populations. Safety outcomes were analysed in the intention-to-treat population. This trial was registered with the Netherlands Trial Register, NL5946. Findings: Between April 12, 2017, and June 3, 2021, 13 267 patients were screened and 1066 were randomly assigned, 533 to each group. 31 were excluded from intention-to-treat analysis of the 2-day group and 30 from the 5-day group owing to errors in recruitment or consent. Appendicectomy was done laparoscopically in 955 (95%) of 1005 patients. The telephone follow-up was completed in 664 (66%) of 1005 patients. The primary endpoint occurred in 51 (10%) of 502 patients analysed in the 2-day group and 41 (8%) of 503 patients analysed in the 5-day group (adjusted absolute risk difference 2·0%, 95% CI −1·6 to 5·6). Rates of complications and re-interventions were similar between trial groups. Fewer patients had adverse effects of antibiotics in the 2-day group (45 [9%] of 502 patients) than in the 5-day group (112 [22%] of 503 patients; odds ratio [OR] 0·344, 95% CI 0·237 to 0·498). Re-admission to hospital was more frequent in the 2-day group (58 [12%] of 502 patients) than in the 5-day group (29 [6%] of 503 patients; OR 2·135, 1·342 to 3·396). There were no treatment-related deaths. Interpretation: 2 days of postoperative intravenous antibiotics for complex appendicitis is non-inferior to 5 days in terms of infectious complications and mortality within 90 days, based on a non-inferiority margin of 7·5%. These findings apply to laparoscopic appendicectomy conducted in a well resourced health-care setting. Adopting this strategy will reduce adverse effects of antibiotics and length of hospital stay. Funding: The Netherlands Organization for Health Research and Development.</p

Small coronary calcifications are not detectable by 64-slice contrast enhanced computed tomography

Recently, small calcifications have been associated with unstable plaques. Plaque calcifications are both in intravascular ultrasound (IVUS) and multi-slice computed tomography (MSCT) easily recognized. However, smaller calcifications might be missed on MSCT due to its lower resolution. Because it is unknown to which extent calcifications can be detected with MSCT, we compared calcification detection on contrast enhanced MSCT with IVUS. The coronary arteries of patients with myocardial infarction or unstable angina were imaged by 64-slice MSCT angiography and IVUS. The IVUS and MSCT images were registered and the arteries were inspected on the presence of calcifications on both modalities independently. We measured the length and the maximum circumferential angle of each calcification on IVUS. In 31 arteries, we found 99 calcifications on IVUS, of which only 47 were also detected on MSCT. The calcifications missed on MSCT (n = 52) were significantly smaller in angle (27° ± 16° vs. 59° ± 31°) and length (1.4 ± 0.8 vs. 3.7 ± 2.2 mm) than those detected on MSCT. Calcifications could only be detected reliably on MSCT if they were larger than 2.1 mm in length or 36° in angle. Half of the calcifications seen on the IVUS images cannot be detected on contrast enhanced 64-slice MSCT angiography images because of their size. The limited resolution of MSCT is the main reason for missing small calcifications

Endothelial shear stress 5 years after implantation of a coronary bioresorbable scaffold

Aims: As a sine qua non for arterial wall physiology, local hemodynamic forces such as endothelial shear stress (ESS) may influence long-term vessel changes as bioabsorbable scaffolds dissolve. The aim of this study was to perform serial computational fluid dynamic (CFD) simulations to examine immediate and long-term haemodynamic and vascular changes following bioresorbable scaffold placement. Methods and results: Coronary arterial models with long-term serial assessment (baseline and 5 years) were reconstructed through fusion of intravascular optical coherence tomography and angiography. Pulsatile non-Newtonian CFD simulations were performed to calculate the ESS and relative blood viscosity. Time-averaged, systolic, and diastolic results were compared between follow-ups. Seven patients (seven lesions) were included in this analysis. A marked heterogeneity in ESS and localised regions of high blood viscosity were observed post-implantation. Percent vessel area exposed to low averaged ESS (7 Pa) did not significantly change (moderate ESS: 76.93% vs. 80.7%, P = 0.546; high ESS: 7.15% vs. 14.31%, P = 0.281), leading to higher ESS at follow-up. A positive correlation was observed between baseline ESS and change in lumen area at 5 years (P < 0.0001). Maximum blood viscosity significantly decreased over 5 years (4.30 ± 1.54 vs. 3.21± 0.57, P = 0.028). Conclusion: Immediately after scaffold implantation, coronary arteries demonstrate an alternans of extremely low and high ESS values and localized areas of high blood viscosity. These initial local haemodynamic disturbances may trigger fibrin deposition and thrombosis. Also, low ESS can promote neointimal hyperplasia, but may also contribute to appropriate scaffold healing with normalisation of ESS and reduction in peak blood viscosity by 5 years